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We propose a unified approach for describing the carrier localization, metal-insulator transitions 
(MITs) and stripe formation in inhomogeneous hole-doped cuprates. The ground-state energy of 
a carrier interacting with a defect and with lattice vibrations is calculated variationally within 
the continuum model and adiabatic approximation. At low doping levels, hole carriers in La-based 
systems with large-radius dopants are localized near the dopants with the formation of hydrogen-like 
impurity centers. As the doping increases, the carriers are liberated from the hydrogen-like impurity 
centers and are self-trapped in a defect-free deformable lattice with the formation of intrinsic large 
polarons. In La-based cuprates with small-radius dopants, hole carriers are self-trapped near the 
dopants with the formation of non-hydrogen-like impurity centers or extrinsic large polarons. The 
charge ordering and formation of different superlattices and energy bands of dopants and intrinsic 
large polarons at their inhomogeneous spatial distribution trigger the MITs in cuprates. We analyse 
the validity of the criteria for the Mott and Anderson MITs in La-based cuprates with hydrogen- 
like impurity centers and show that such MITs in these systems are unlikely possible. Using the 
uncertainty principle we derive the quantitative criteria for the new MITs caused by strong carrier- 
dopant-phonon and carrier-phonon interactions. We show that the MITs in La-based cuprates with 
large-radius dopants are driven by intrinsic self-trapping and charge ordering, while the anisotropy 
of dielectric constants fosters the MITs accompanied by the formation of static stripes (in carrier- 
poor domains) and dynamic ones (in carrier-rich domains) in the range of doping x ~0. 04-0. 125. 
The so-called 1/8 anomaly is especially pronounced in these systems. The small-radius dopants 
and the anisotropy of dielectric constants favor carrier localization, MITs and stripe formation in 
a wide range of doping (including also "magic" doping x=l/8) in other La-based cuprates, where 
x=l/8 is no longer "magic" doping. Finally, the MIT and 1/8 anomaly in cuprates caused by 
the commensurate ordering of planar large bipolarons and their condensation into a liquid are also 
discussed. The obtained results are in quantitative agrement with the experimental data on polaron 
formation, MITs and stripe formation in La-based cuprates. 

PACS numbers: 71.30.+h, 71.38.-H, 71.55.-i, 74.20.Mn, 74.72.-h, 74.72.Dn 
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I. INTRODUCTION 



The carrier localization, metal-insulator transitions 
(MITs) and stripe formation in hole-doped cuprates 
are the most challenging problems in condensed matter 
physics. These phenomena have attracted much atten- 
tion because they are closely related to high-T c supercon- 
ductivity in cuprates 1 -"— . The mechanisms for carrier lo- 
calization and MITs in the cuprates and related materials 
have been under discussion since the discovery of high- 
Tc superconductivity- ^ 10 ' 12 ! 15 ' 17- — . Another important 
problem of the physics of high-T c cuprates is the role of 
the electronic inhomogeneity and charge ordering in the 
stripe formatio n 11 ' 14 ' 20 ' 21 , which is intimately related to 
carrier localization and MITs in these materials. For 
cuprates within the insulating antiferromagnetic (AF) 
phase, the electron correlation is likely to dominate the 
carrier localization; when a MIT is approached, it is not 
obvious which process will dominat o 12 ' 22 . Actually, the 
MITs in the cuprates might be different from the Mott 
and Anderson transitions. In these systems, the phenom- 
ena of carrier localization seem to be much more com- 
plicated by electron-defect-lattice and/or electron-lattice 
interaction effects which become very pronounced close 
to the MIT— 122. Generally, the electron-defect-lattice 



and/or electron-lattice interactions and the electronic in- 
homogeneity are known to foster carrier localization in 
solids. In particular, the strong electron-defect-phonon 
and electron-phonon interactions foster carrier localiza- 
tion by causing extrinsic and intrinsic self-trapping. At 
the same time the gross electronic inhomogeneities might 
play an important role in carrier localization, MITs 
and formation of insulating and metallic stripes in the 
cuprates. There is now increasing experimental evidence 
for the existence of extrinsic polarons (i.e., charge carri- 
ers self-trapped near the dopants or impurities), intrin- 
sic polarons (or charge carriers self-trapped in a defect- 
free deformable lattice) and electronic inhomogeneities 
in these materials ^ 20 ' 21 ' 23 ' 24 . In addition, various ex- 
periments on hole-doped La-based cuprates show that 
the MITs and stripe formation occur in underdoped (at 
doping levels x — 0.05 — 0.125)^^, optimally doped 
(x = 0.15-0.16) 8 , and overdoped (x = 0.16-0.25) 16 i 21 i 27 
regimes. However, the role of electronic inhomogeneities 
and strong carrier-defect-phonon and carrier-phonon in- 
teractions, which is known to be very important in the 
vicinity of the MIT a 20 ' 22 ' 28 , has not been fully explored 
in the underdoped regime. While the effects of ionic (or 
dopant) size and anisotropy of dielectric constants, which 
may play a key role in carrier localization and MITs, 
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have not been studied theoretically in inhomogeneous 
La-based cuprates in a wide hole concentration range 
from x ~ 0.02 (lightly doped regime) to x ~ 0.25 (heav- 
ily overdoped regime). In this work we examine such 
inhomogeneous systems in which the electronic charge 
inhomogeneities in conjunction with the strong carrier- 
defect-phonon and/or carrier- phonon interactions induce 
the carrier localization and MITs accompanied by the 
stripe formation in a wide range of doping. In so do- 
ing, we study the MITs and stripe formation in high- 
T c cuprates driven by the inhomogeneous local doping, 
charge ordering and self-trapping of hole-carriers near de- 
fects (impurities) or in a defect-free deformable lattice. 

This paper is organized as follows. First, in Sec. II, we 
investigate the extrinsic (or defect-assisted) and intrinsic 
self-trapping of carriers, and analyze the possibility of 
formation of extrinsic (hydrogen-like and non-hydrogen- 
like) and intrinsic self-trapped states in hole-doped La- 
based cuprates. Then, in Sec. Ill, we consider the possi- 
ble forms of charge ordering and the formation of differ- 
ent superlattices and energy bands of dopants (i.e., accep- 
tors) and intrinsic large polarons at their inhomogeneous 
spatial distribution. In Sec. IV, we examine the validity 
of the criteria for the Mott and Anderson transitions in 
La-based cuprates with hydrogenic impurity centers. We 
argue that such transitions in the cuprates are unlikely 
possible. We give our view of the new MITs in La-based 
cuprates with extrinsic large polarons (non-hydrogen-like 
impurity centers) or intrinsic large polarons and derive 
the quantitative criteria for the MITs by using the perti- 
nent uncertainty relation. We show that in these systems 
the extrinsic and/or intrinsic self-trapping of hole carriers 
and the distinctly different charge ordering result in the 
new MITs which occur at different doping levels. Fur- 
ther, we find that, the ionic radius of dopants and the 
anisotropy of dielectric constants strongly influence on 
carrier localization, MITs, stripe formation and so-called 
1/8 anomaly in La-based cuprates. These predictions are 
consistent with distinctive features of MITs, stripe for- 
mation and 1/8 anomaly observed in the cuprates. Fi- 
nally, we briefly discuss the other type of MIT and the 
1/8 anomaly in hole-doped cuprates caused by the com- 
mensurate ordering of planar large bipolarons and their 
condensation into a liquid. The paper concludes with 
Sec. V, in which the principle results are summarized. 



II. EXTRINSIC AND INTRINSIC 
SELF-TRAPPING OF HOLE CARRIERS 



assumed to be within both a three-dimensional (3D) and 
two-dimensional (2D) deformable medium, the last one 
being CuC-2 layers 1 ^. The large static dielectric constant 
Eq (>> 10) both parallel and perpendicular to the CuC>2 
layers 1 ^ indicate that lattice ions may be significantly 
displaced in all three directions. In reality, the CuC>2 
based materials may be approximated as a 3D deformable 
medium 1 -. As it was pointed out in Ref^ -, the existence 
of the 3D bismuth oxide high-T c superconductors, which 
belong to the same class of materials as the copper ox- 
ides, casts an extra doubt on the theories based mainly 
on the low dimensionality of the crystal structure. There 
is also convincing experimental evidence that the consid- 
eration of cuprates as 3D systems may appear to be more 
appropriate (see Ref 12 ' 31 ). In these anisotropic 3D polar 
materials, the hole carriers interacting both with lattice 
vibrations (i.e., acoustic and optical phonons) and with 
lattice defects (e.g., dopants or impurities), can easily 
be self-trapped near the lattice defects. Such extrinsic 
self-trapping of hole carriers leads to the formation of 
a moderately deep extrinsic polaronic or a shallow hy- 
drogenic state in the CT gap of the parent cuprates. 
The theory of carrier self-trapping has been developed 
in Refs^£ & 32 ' 33 in the framework of the continuum ap- 
proach. As shown in Refs i 33 ' 34 , the continuum theory of 
carrier self-trapping in a defect-free system can be eas- 
ily extended to the system consisting of a defect and a 
bound charge carrier. Such a continuum approach is bet- 
ter suited for the quantitative analysis. In the framework 
of the continuum model, we make an attempt to analyze 
the possibility of the formation of hydrogen-like and non- 
hydrogen-like impurity states as well as polaronic ones in 
the CT gap of the cuprates. 

The ground-state energy of a hole carrier in 3D polar 
crystals can be calculated variationally in the continuum 
model and adiabatic approximation, taking into account 
the short- and long-range carrier-defect-phonon interac- 
tions. The total energy of the interacting carrier-defect- 
phonon system is given by the functional 34 
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The undoped layered cuprates are typical charge- 
transfer (CT)-type Mott insulators. According to the 
band calculations and spectroscopic dat a 10 ' 29 , the elec- 
tronic structure of these parent CuCvbased compounds 
is well described by the three-band Hubbard model and 
the oxygen valence band lies within the Mott-Hubbard 
gap. Upon hole doping, the oxygen valence band of the 
cuprates is occupied by free holes. The hole carriers are 



where ip( r ) an d rn* are the wave function and effective 
mass of a carrier, e — £oo/(l — v) is the effective dielectric 
constant, rj = £oo/£o> £ oo is the high frequency dielectric 
constant, K is an elastic constant, Ed and EdD are the 
deformation potentials of the carrier and the defect, re- 
spectively, Vq is the short-range defect potential, Z is the 
charge state of the defect. 

Now, we may choose the trial wave function ip(r) in 



3 



the Gaussian form 



%jj{r) = Nexp [— {crrf 



(2) 



and make a variational calculation of Eq. (JT]) with re- 
spect to the parameters contained in ip( r ), where N — 
(2/7r) 3 / 4 a 3 / 2 is the normalization factor, a — y^fi/aa is 
the variational parameter, /3 — ao/a is the localization 
parameter, a is the radius of localization limited between 
ao (the lattice constant) and oo (free carrier). Inserting 
Eq. @ into Eq. (|T|) and performing the r integration in 
Eq. ([I]), one can obtain the following functional 



E{p) =B[/3 2 - 3l /3 3 - 52 /3] 



(3) 



where B = 3ttH 2 /2m*a 2 ., g 1 = (E%/2Ka%B)(l + b s ), 
b s = 2 5 / 2 [E dD /E d -KV /Ej], 52 = (e 2 / £oo a B)(l - 
V + b l r 1 ),b l = 2 3 / 2 Z. 

When 3<7i<72 < 1 the functional (j3|) has a minimum 
at (3 m in = [l — \/l — Sg\g2\ /3gi which corresponds to 
the formation of a large-radius impurity state or extrin- 
sic polaronic state in the CT gap of the doped cuprates. 
A carrier interacting with lattice vibrations becomes a 
large polaron with a polarization (or deformation) cloud 
extending over a wide region. We can determine the 
ground-state energy of such an intrinsic large polaron 
from Eq.© at b s =0 and Z = 0. At 3gig 2 < 1 the po- 
tential barrier height separating large- and small-radius 
self-trapped states of a hole carrier is equal to 34 i 35 
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While the ground-state energy of a self-trapped carrier 
is given by 3 ^ 
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The signs of Ed and E d D for holes and small-radius 
defects always positive, while E d ]j for large-radius de- 
fects is negative 3 ^. In La2- x Sr x Cu04 (LSCO) the ra- 
dius of Sr 2+ ions is larger than that of La 3+ ions2£, so 
that for SV 2+ ion Z = 1, E dD < and b s < 0. In 
this case, the short-range part of the impurity potential 
in Eq. ^ is repulsive and the hole-lattice interactions 
near such dopants are suppressed by this repulsive de- 
fect potential. At a weak carrier-lattice interaction the 
localized impurity state may have hydrogen-like charac- 
ter described by a rigid lattice model^ 8 .. Therefore, in 
order to study the Mott MIT in LSCO, we can con- 
sider the hydrogenic impurity centers having the Bohr 
radius an — 0.529eo(m e /m*)A and the ionization energy 
Ej = e 2 /2eoa,H in the Hubbard model, where m e is the 
free electron mass. However, the substitution of small- 
radius cations (e.g., Ca 2+ and Nd 3+ ions) for La 3+ ions in 
La2Cu04 and for Sr 2+ ions in LSCO leads to the forma- 
tion of the non-hydrogenic impurity centers or extrinsic 
large polarons for which the Hubbard model is already 
inapplicable. The distinctive feature of the cuprates is 



their very large ratio of static to high frequency dielec- 
tric constants 1 "! 3 ^. This situation is favorable for carri- 
ers attracted to polarization well created by the other 
ones or to Coulomb centers (dopants) to form 3D intrin- 
sic or extrinsic large bipolarons. Such bipolarons with 
small binding energies can be formed in the lightly doped 
cuprates at 77 < 0.1— and they become unstable in the 
underdoped regime. 



III. POSSIBLE TYPES OF CHARGE 
ORDERING AND FORMATION OF DIFFERENT 
SUPERLATTICES 

The hole-doped cuprates are inhomogeneous or disor- 
dered systems, where the dopants and charge carriers are 
distributed inhomogeneously. Such dopant-driven and 
carrier-driven inhomogeneities may produce regions with 
different doping levels. The electronic inhomogeneity 
commonly exists in high-T c cuprates regardless of doping 
level and the underdoped cuprates become more inhomo- 
geneous than overdoped ones 2 ^ Inhomogeneous distri- 
butions of carriers result from their interaction with one 
another (e.g. phase separation) and/or from their in- 
teraction with an inhomogeneous distribution of dopants 
or micro-structural defects such as dislocations 39 . One 
can assume that the electronic disorder or inhomogene- 
ity in the cuprates leads to the charge segregation into 
carrier-rich and carrier-poor regions. As the doping level 
increases towards underdoped region, specific charge or- 
dering takes place in these regions and distinctly differ- 
ent superlattices and energy bands of dopants (impurities 
with trapped free carriers or large polarons) and self- 
trapped carriers (intrinsic large polarons) are formed at 
their inhomogeneous spatial distribution. In particular, 
the hydrogenic impurity centers (impurities with loosely 
bound free carriers or large polarons) and non-hydrogenic 
ones (impurities with tightly bound large polarons) are 
assumed to form the superlattices with the lattice con- 
stant ai and coordination number z. The charge order- 
ing in carrier-poor and carrier-rich domains results in the 
formation of simple cubic, body-centered cubic and face- 
centered cubic superlattices with coordination numbers 
z =6, 8 and 12, respectively. In this case the formation 
of different impurity bands in the cuprates is described 
by the tight-binding approximation and the widths of the 
impurity bands can be determined from the relation 



Wi = 2zJi, 



(6) 



where Jj — h 2 /2m*joj is the hopping integral between 
nearest-neighbour impurity centers, m*j is the effective 
mass of charge carriers in the impurity band. 

Further, we believe that the intrinsic large polarons 
just like impurity centers form different superlattices with 
the lattice constant a p and the widths of the polaronic 
bands just as the widths of the small-polaron bands^S are 
determined from the expression 

W p = 2zJ p , (7) 
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where J p = h 2 /2m*a 2 is the hopping integral between 
nearest-neighbour sites of the polaronic superlattice, m* 
is the effective mass of large polarons. 

At low doping level the extrinsic (or intrinsic) large 
polarons form a simple cubic superlattice with aj or 
a p > 2R P , where R p is the radius of the polaron. The 
width of the impurity band decreases with increasing 
inter-dopant separation, so that the impurity-band ef- 
fective mass mi increases with inter-dopant separation or 
with decreasing doping level. Clearly, in the lightly doped 
regime the narrow impurity and/or polaronic bands are 
formed in the CT gap of the parent cuprates and the en- 
ergy gaps exist between the oxygen valence band and the 
impurity or polaronic bands. These energy gaps (or pseu- 
dogaps) disappear in sequence when the polaronic effects 
disappear in carrier-rich and carrier-poor domains in the 
overdoped regime. 

We now make some remarks about the confinement en- 
ergy of a large polaron and the width of the large-polaron 
band in the cuprates. The large-polaron confinement en- 
ergy is about h 2 /rripR p 39 ' 41 . When the Frohlich electron- 
phonon coupling constant a is large (a ^> 1) the effective 
mass of a large polaron is given by the formula 4 ^ 

16to* 4 

m P = (8) 

The Frohlich coupling constant a for LSCO is esti- 
mated to be a = 5. 7 39 i 42 . For m* = mj^, we then ob- 
tain from (8), m p ~ 2.14m e , which is close to the mea- 
sured value of m p — 2m e in LSCO 9 . Using the Bohr 
radius of the polaron R p = 0.529eo(m e /m p )A ~ 7 A2A 
(at Eq = 30) we obtain H 2 /m p R 2 ~ 0.065eU, which is less 
than E p (see Sec. IV) and is larger than the characteristic 
phonon energies hu) = 0.04— 0.06eU in the cuprates ^ 23 ' 43 . 
As the doping level is increased to underdoped levels 
at which the insulator-to-metal transition takes place, 
the distance between large polarons decreases and they 
can form closely packed simple cubic superlattice with 
a p = 2R p . In this case the large-polaron bandwidth is 
about W p ~ O.leV (i.e., W p > E p ). In contrast, the 
small-polaron band will be very narrow (W p <C froS) since 
the effective mass of a small polaron is much larger than 
a free-electron mass. 



IV. METAL-INSULATOR TRANSITIONS AND 
STRIPE FORMATION 

Now, we discuss how the extrinsic and/or intrin- 
sic self-trapping of carriers, the dopant- and carrier- 
driven charge inhomogeneities and the charge ordering 
are related to the MITs and stripe formation in hole- 
doped cuprates. Various experiments have confirmed 
that for light doping, La2CuC>4 behaves like a conven- 
tional semiconductor 9 . As described above, large-radius 
dopants or impurities with E^d < and b s < may form 
the hydrogenic acceptor centers in La-based cuprates 
and the Hubbard model based on the strong on-center 



Coulomb repulsion U is applicable for studying the MIT 
in these systems. Then, the Mott MIT point is deter- 
mined from the condition Wi — Ej — 1.16J7 which was 
used in the derivation of the following criteria for the 

Mott MITs 19 : n c /3 a H ~ 0.289 (z = 6), nJ 3 a H ^ 0.315 

1/3 

(z = 8) and n c au — 0.324 (z — 12), where n c is the 
critical concentration of the hydrogenic impurities initiat- 
ing the Mott MIT in doped systems. If we take m* — m e 
(for free carriers) and Eq = 30, we find an = 15.87A in 
LSCO. According to the above Mott criteria, the MITs 
would occur at the dimensionless hole concentrations 
x c = n c /n a ~ 0.00115 - 0.00161, where n a = l/V a is 
the density of the host lattice atoms, V a ~ 190A 3 is 
the volume per Cu02 unit in the orthorhombic LSCO. 
The above values of much smaller than the value 

of x c ~ 0.02 at which the destruction of the antiferro- 
magnetic (AF) order is observed in LSCO^. However, 
the large-radius dopant and large polaron may form the 
hydrogenic impurity center in LSCO. When b s < (or 
EdD < 0), the defect and the self-trapped hole repel each 
other at short distance but attract each other at long 
distance. Therefore, the hydrogenic self-trapped state of 
a large polaron in LSCO has an effective Bohr radius 
an = 0.529£o(m e /m p )A The existence of such hydro- 
genic acceptor centers in La-based cuprates is now exper- 
imentally well established 9 . 

According to Mott and Davis 3 ^, the Mott's criterion 
for the MIT can be used as the condition for the metallic 
behavior of a degenerate polaron gas if the Bohr radius 
an of polarons is larger than the interatomic distance. If 
we take m p = 2m e and So = 27, we find an — 7.14A, 
so that the above Mott criteria for the MITs can be 
applied to the cuprates. In this case the Mott transi- 
tions would occur at x c —0. 0112-0. 0178. The quantita- 

1 /3 

tive criteria for the Anderson MITs, n c an — 0.289 (for 
z=6), n c /3 a H ~ 0.293 (for z=8), and n c /3 a H ~ 0.272 
(for z=12), derived in Refi^ 9 - predict nearly the same val- 
ues of x c . We see that the Mott transition (for shal- 
low impurity centers) would occur at a much smaller 
dopant concentration than the critical doping concen- 
tration for the MIT in the cuprates. For example, the 
MITs in LSCO are observed at more higher doping levels 
x c ~ 0.05 — O.O?-^^^^ and these experimental results 
cannot be reconciled with the above criteria for the Mott 
MITs. In our opinion, the hole carriers liberated from 
the large-radius dopants or hydrogenic impurity centers 
a t %c <^ 0.02 are constrained to remain away from these 
impurities and they are self-trapped in a defect-free de- 
formable lattice with the formation of strong-coupling 
intrinsic large polarons, which can acquire itineracy at 
higher doping levels 2 ^. 

Another possibility is that small-radius dopants (e.g., 
Ca 2 + and Nd 3+ ) with E dD > (or b s > 0) favor the for- 
mation of the extrinsic large polarons (non-hydrogenic 
impurity centers) in La-based cuprates. In the present 
case, both the short and long range parts of the carrier- 
defect interaction is attractive, so that hole carriers are 
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self-trapped near such dopants. Polaronic effects will be 
stronger near the small-radius dopants, leading to car- 
rier localization over a broader range of doping. We can 
evaluate the possibility for the existence of a potential 
barrier separating large- and small-radius intrinsic pola- 
ronic states in 3D hole-doped La-based cuprates using 
the relation Q. 

According to the spectroscopic data, the Fermi energy 
of the undoped cuprates is about Ep ~ 7e V 44 ^ 45 . To de- 
termine the value of the short-range carrier-phonon cou- 
pling constant g\, we estimate Ed as Ed = (2/3)Ep. The 
values of other parameters are K = 1.4 • 10 12 dyn/cm 2 — , 
m* = mJL, £oo = 3 — 5i, ao — QA (for orthorhombic 
LSCO), b s = and Z = 0. Then B = 1 eV, g 1 ~ 0.0576 
and g2 — 2.4(1 — rfi/soQ. For the cuprates, typical values 
of i] range from 0.02 to 0.122. Using Eq.Q, we obtain 
E a value of ~ 38 eV at = 4 and rj = 0.06. As al- 
ready seen, the large- and small-radius polaronic states 
are separated by very high potential barrier. This bar- 
rier prevents the formation of small polarons in the 3D 
cuprates. At b s = 0.5 and Z = 1 the height of such a 
potential barrier separating large- and small-radius ex- 
trinsic polaronic states is ~ 14 eV. It follows that the 
relevant charge carriers in the anisotropic 3D cuprates 
are the large intrinsic and extrinsic polarons (or bipo- 
larons). While the small polarons and bipolarons may 
be formed in the CuC>2 layers of the cuprates, where the 
self-trapping of carriers may occur more easily due to the 
absence of the potential barrier between the large- and 
small-radius self-trapped states in 2D system s 47 ! 48 . How- 
ever, such (bi)polarons having very large effective masses 
and narrow energy bands tend to be localized rather than 
mobile. 

Small or large polarons may be expected in polar ma- 
terials with a small or large hopping integral J between 
nearest-neighbour sites of the crystal lattice. The ra- 
dius of the polaron R p is of order ao(2J/E p ) and the 
formation of a polaron whose radius is large compared 
to the lattice constant do requires that 2J/E p 1 (see 
also Ref42) or W p /E p > In hole-doped cuprates 

(LSCO, YBa 2 Cu 3 07-5 and Bi 2 Sr 2 CaCu208+5) large po- 
larons have the effective masses m p ~ (2 4)m<r^i. 
Therefore, Wj and W p may be comparable with the bind- 
ing energies E p and E p of the extrinsic and intrinsic large 
polarons since m p , mj(> m p ), aj and a p would decrease 
with increasing carrier concentration. At m/ < 2m e 
and m p < 2m e the magnitudes of Wi and W p may be- 
come larger than E l v and E p , which are determined from 
Eq.© by using the above presented values of parame- 
ters. The obtained results are summarized in Table 1. 
When 7] increases, the ionization energy E p of the non- 
hydrogenic impurity centers (with Z = 1) increases. On 
the contrary, the binding energy of the intrinsic large po- 
larons E p decreases. Interestingly, E p and E p increases 
markedly with decreasing from 5 to 3. We see that 
the value of E^—OAS eV obtained at 77=0.10 and Z=l is 
consistent with the experimental data for lightly doped 
La2Cu04+j^ and the impurity band observed in these 



TABLE I: Binding energies of intrinsic and extrinsic large 
polarons for Soo = 4, b s = 0.5 and different values of 77. 



n 


0.00 


0.02 


0.04 


0.06 


0.08 


0.10 


0.12 


E p , eV 


0.0916 


0.0879 


0.0844 


0.0809 


0.0774 


0.0741 


0.0708 


z = 


E J P , eV 


0.0925 


0.0888 


0.0851 


0.0816 


0.0781 


0.0747 


0.0714 


Z = 1 


E l v , eV 


0.0925 


0.0995 


0.1067 


0.1142 


0.1220 


0.1301 


0.1384 



systems at 0.13 eV is associated with the extrinsic large 
polarons rather than the small polarons bound to impu- 
rities. 

As we have discussed above, the extrinsic large po- 
larons bound to dopants may also form different super- 
lattices in inhomogeneous hole-doped cuprates and the 
energy bands (which are formed in the lightly doped 
cuprates) of such polarons may exist, thus permitting 
charge transport by means of intra-band conduction. In 
case of narrow polaronic bands, charge transport be- 
comes hopping-like and is caused by intra-band hopping 
processes^. One can assume that if the bandwidth of 
the extrinsic large polarons exceeds some critical value, 
their intra-band conduction becomes metal-like. We at- 
tempt to find the new criteria for such MITs in doped 
materials. 

When the carrier concentration is low, the carriers oc- 
cupy low- lying bound states (i.e., extrinsic or intrinsic 
polaronic states) first. As soon as all the localized states 
were filled at some critical carrier concentration n = n c 
(or x — x c ), the carriers start to occupy the itinerant 
states and the transition of the system from the insulat- 
ing phase to the metallic one occurs. The possible values 
of n c at which extrinsic or intrinsic large polarons forming 
different superlattices acquire itineracy are determined 
from the appropriate criteria for the MITs. The condi- 
tions for carrier localization or delocalization can be ob- 
tained by using the uncertainty principle: ApAx > h/2, 
where Ap and Ax are the uncertainties in the momentum 
and coordinate, respectively. This uncertainty relation 
can be written as 



Ax-AE 



h 2 {Ak) 2 
2m* 



1 



2Afc' 



(9) 



where AE and Ak are the uncertainties in the energy and 
wave vector, respectively. The expression h 2 (Ak) 2 /2m* 
in Eq. ([9]) represents the uncertainty in the energy of 
free carriers. Taking into account that the uncertainties 
in the energy and wave vector in the impurity band are 
about Wj and 1/oj, respectively, the relation ^ can 
be rewritten in the form (cf. a similar relation formerly 
derived by Ridley 4 ^ using the uncertainty principle) 

Ax ■ AE =■ Wiai/2 (10) 

On the other hand, the uncertainty in the energy AE 
of the localized state of carriers is of the order of E^, 
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whereas the uncertainty in the coordinate Ax of the ex- 
trinsic large polarons is of the order of the polaron radius 
Ri. Then, the condition for the setting-in of carrier lo- 
calization can be written as 



K R i Z Wiai/2, 



(11) 



so that the new criterion for the MIT can be derived from 
this condition and written as 



Wi 



0.5 



(12) 



which is similar to the results for analogous prob- 
lems of Holstcin localization, Mott and Anderson MITs. 
In particular, localization in the one-dimensional Hol- 
stein molecular-crystal model is governed by the ratio 
of the local polaron binding energy to the inter-site 
bandwidt h 39 ! 53 . Similarly, Anderson's localization in An- 
derson's disorder model is governed by the ratio of the 
local-site energy (or the random potential Vq) to the elec- 
tron bandwidth and Mott's MIT is governed by the ratio 
of the Hubbard on-site Coulomb interaction energy U 
to the V F 38 i 41 ' 52 . For intrinsic large polarons, instead of 
(p~2|) . we could write 



E p 



0.5 



Rn 



(13) 



With increasing the inter-polaron distance R, the po- 
laron band is continuously narrowed, finally ending in 
the discrete levels (at R 3> R p ) of the non-interacting 
polarons. On the contrary, with decreasing R, the large 
polarons begin to interact one another through the over- 
lap of their structural distortions and the polaronic states 
are broadened into an energy band. As long as the dis- 
tance between nearest-neighbour polarons is larger than 
2R p and the overlap of the deformation clouds of polarons 
is small, one would expect the large-polaron bands to be 
very narrow. In this case large polarons are localized or 
confined to their potential wells and the system is con- 
verted into an insulator. For R = 2R P (which is assumed 
to correspond to the packing of spheres with radii equal 
to the polaron radius), the deformation clouds of large 
polarons begin to overlap strongly and the large polarons 
are delocalized and characterized by metal-like transport 
in sufficiently broadened energy bands. 

If the extrinsic large polarons form such closely packed 
simple cubic, body-centered cubic and face-centered cu- 
bic superlattices with aj = 2Rj (z = 6), a/ = (4/\/3)-Rr 
(z = 8) and ai = (2\/2)Rj (z — 12), the appropriate den- 
sities of such carriers per unit cells of the superlattices 
are n = 1/a 3 , n = 2/a 3 and n = 4/a 3 , respectively. For 
these cases, three successive MITs are feasible and the 
localized extrinsic large polarons in the cuprates start to 
occupy the itinerant states at different hole concentra- 
tions determined from relation (fT2|) . This means that 
the system gradually changes from the insulating phase 
to the metallic one. When the extrinsic large polarons 
form simple cubic and face-centered cubic superlattices, 



such insulator-to-metal transitions occur at Wj = E d and 
Wi = Ep/y/2, respectively. We argue that the electronic 
inhomogeneities and the ordering of polaronic carriers 
with the formation of different superlattices can produce 
various types of self-organized electronic structures in the 
form of stripes. Using the relation (|T2"j) . we obtain the 
following criteria for the new MITs: 



zh 2 



.3/2 



for z = 6, 



71 ^ zh 2 



3/2 



for 



.5/4 f ^Ef 



3/2 



for 



12. 



(14) 



(15) 



(16) 



We are now in a position to evaluate n c in La-based 
cuprates using Eqs. (fT4)l . (fT5)) and (fT6)l . The effective 
masses of carriers in cuprates deduced from ARPES and 
electronic specific heat data (see Reb 5 ^) at different dop- 
ing levels are slightly different and equal to 2.1 - 2.5 times 
the free electron mass. Therefore, we can evaluate n c 
by taking mj — 2.5m e for La2_ x Ca x Cu04. Then for 
f] = 0.02 and E* = O.lOeV we find x c ~ 0.065 - 0.080. 
When £oo = 4, rj = 0.02 - 0.12 and E* = 0.100 - 0.138 
eV (see Table 1), the MITs and stripe formation occur 
in these systems at x c ~ 0.065 — 0.131. At = 3 
and rj = 0.06 we obtain E p ~ 0.2eV. Then the metal- 
insulator boundary of La2- x Ca x Cu04 lies in the over- 
doped regime at x c ~ 0.226 (for z = 8). 

Further, the double substitution of smaller cations 
for host lattice ions and dopants may also favor the 
MITs and stripe formation occurring in a wide range 
of doping of La-based cuprates. When the doping in- 
creases, there is a significant probability of at least two 
neighbouring La 2+ ions replaced by large-radius Sr 2+ 
ion (with Z = 1, V Sr > and E^ < 0) and by 
small-radius Nd 3+ ion (with Z = 0, V Nd > and 
E dD > °) in La 2 - x _ y Nd y Sr x Cu0 4 (where y » x). 
The attractive potentials of these two-dopant centers 
may be greater than those of more separated dopants. 
For La2- x - y Nd y Sr x Cu04 with such two-dopant centers, 
we should replace Vb and EdD by V$ r + V^ d and 
EdB + EdS m Eq. JTJ) or in the expression for b s . At 
E dD + E dD > ( Le -> b s > 0) th e short-range part of 
the two-dopant potential is attractive and the hole car- 
riers are self-trapped near the two neighbouring dopants 
with the formation of the extrinsic large polarons. The 
two-dopant-driven charge inhomogeneity and ordering 
lead to formation of the superlattices of such extrin- 
sic large polarons. In order to illustrate the effects 
of anisotropy of Soa and the short-range defect poten- 
tial on x c , we show in Fig.l results of our calculations 
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FIG. 1: Critical carrier concentration x c as a function of Ex at 
z = S, rj = 0.06, b s = 0.1 and b s = 1.5. 



FIG. 2: Critical carrier concentration x c as a function of r) at 
z = 8, b s = 0.5 and £oo = 4. 



for hole-doped La-based cuprates with two-dopant cen- 
ters for which Z = 1 and b s > 0. Fig.l shows that 
the anisotropy of Sqo and the size effect of small-radius 
dopants are the main driving forces for the new MITs 
and stripe formation in a wide range of doping (includ- 
ing also x = 1/8) in La-based cuprates containing two 
types of dopants (i.e., large- and small-radius dopants). 
Indeed, in hole-doped cuprates La2- x - y Nd y Sr x Cu04 and 
La2_x-yEu y Sr x Cu04 with small-radius dopants Nd 3+ 
and Eu 3+ , the static stripe phases are observed in a wide 
range of carrier concentration and not restricted to a nar- 
row range around x — 1/8 in which superconductivity is 
suppresse d 16 ^ 7 ' 55 . 

As mentioned above, the size of impurities is also a 
key parameter in the cuprates containing large-radius 
dopants and the formation of intrinsic large polarons in 
these systems becomes possible at distances far from the 
impurities. Such materials are the compounds LSCO and 
La2- x Ba x Cu04 (LBCO) in which large polarons form 
different superlattices and acquire itineracy at some crit- 
ical carrier concentrations n c determined from the rela- 
tions (fill) . (|T5|) and ([TFj)) by replacing mi and by m p 
and E p in these relations. Using the values of rj = 0.02, 
Eoo = 4, m p = 2.1m e and E p ~ 0.088e\^ (see Table 1) 
for LSCO and LBCO, we find x c ~ 0.049 (for z=6) and 
x c ~ 0.051 (for z=8). These values of x c are in excel- 
lent agreement with the experimental value x c = 0.05 in 
LSCOi^. If we take into account a possible anisotropy 
of Eoo =3 — 4 and rj = 0.02 — 0.10, then we obtain 
x c ~ 0.04 — 0.125, which are also well consistent with ex- 
isting experimental data on MITs and stripe formation in 
LCSO and LBCO (see Refs^^^^). In Fig. 2 we plot 
the dependence of x c on rj for z = 8, b s = 0.5 and Soo = 4. 
Fig. 2 shows that the value of x c increases markedly 
with increasing rj for the cuprates containing large-radius 
dopants. Our results clearly indicate that x = 1/8 is 
truly "magic" or particular doping level for hole-doped 
cuprates LSCO and LBCO. In these systems, the sup- 
pression of superconductivity should be expected for x 



close to 1/8. However, the doping level x — 1/8 is not 
" magic" doping level for the formation of stripes in other 
hole-doped cuprates, where the effects of small-radius 
dopants on MITs and static stripe formation are strong 
enough. The stripe formation driven by the MITs occurs 
in such systems not only at x — 1/8 but also at much 
more higher doping levels. Indeed, experimental studies 
on La2-xBa x _yM y Cu04 (where M = Sr, Ca) show that 
the 1 /8 anomaly is reduced by substituting smaller diva- 
lent cations instead of large-radius Ba 2+ — . In contrast, 
the ordering of intrinsic large polarons is important for 
the MITs and stripe formation in LSCO and LBCO. Dou- 
ble doping experiments on La2- x - y Ba x Th y Cu04 indi- 
cate that superconductivity is suppressed when the hole 
concentration p — x — y (not the individual or total 
dopant concentrations x, y and x + y) is 1/ 8 56 ' 57 . In 
this system, the Th 4+ ion acts to compensate the neg- 
ative effective charge of Ba-site, so that the hole con- 
centration is now given by p — x — y. This means 
that the presumed superlattice and stripe formation ac- 
companied by the suppression of superconductivity in 
La2- x - y Ba x ThyCu04 are associated with the ordering 
of intrinsic large polarons 39 rather than dopants. 

The formation of 3D and 2D large bipolarons is also 
expected in the cupratesi^i^— . Such bipolarons can 
only exist if the Frohlich coupling constant a is greater 
than the critical value a c and when 77 is smaller than the 
critical value rj c . For instance, for 3D and 2D systems, 
the following critical parameters were obtained by us- 
ing the Gaussian type trial wave functions: a^ D ^ — 5.8, 
■q^ D) ~ 0.13&S2 and a { c D) ~ 2.94, j]f D) ~ 0.167 (for 

a = 5), rii 2D) ~ 0.173 (for a = 10)^. The ratio of 3D 
bipolaron binding energy E^b to twice the polaron en- 
ergy 2E p is about E bB /2E p ~ 0.22 at 77 — > 0§I, whereas 
the ratio of 3D extrinsic bipolaron binding energy E^ B 
to twice the extrinsic polaron energy 2E^ is somewhat 
smaller than the value of E{,b /2E p 62 . Furthermore, the 
binding energies of 3D extrinsic and intrinsic large bipo- 
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larons decrease with the increase of r\. Therefore, such 
bipolarons may exist at low doping levels. As the car- 
rier concentration increases, 3D large bipolarons disso- 
ciate first into two large polarons and then these po- 
larons acquire itineracy at a higher doping level x — x c 
which corresponds to the onset of the metallic behav- 
ior of the system. However, the properties of 2D large 
(bi) polarons may be quite different from those of 3D 
large (bi)polarons. It is clear that the binding energy and 
the effective mass of a bipolaron would increase greatly 
with decreasing dimensionality. In particular, the bind- 
ing energy of 2D large bipolarons E^ 3 ^ is expected to be 

The 2D polaron and bipolaron 



much greater than E^\ 



ground-state energies E^ D> and E^ 1 can be obtained 
approximately using the scaling relation; 



(2D) 



,58.63. 



E^ D \a) 



±E£ D \(3n/4)a) 



and 



E 



(2D) 



(a, U) = -E 



(3D) 



37T 37T 

—a, —U 



(17) 



(18) 



where E^ D ^ and Eg D ^ are the 3D polaron and bipolaron 
ground-state energies, respectively, U is the strength of 
the Coulomb repulsion between the pairing carriers. 

In the strong coupling regime a > a c , the ground-state 
energies of 3D large polaron and bipolaron are propor- 
tional to a 2 —. In this case the 2D large bipolaron binding 
energy is given by 



E. 



(2D) _ 
bB 



l%\ E, 



? (3D) 
J bB 



3.7 E, 



(3D) 
bB 



(19) 



from which it follows that the bipolaron binding energy 
will be much larger in 2D systems as compared to 3D 
ones and, therefore, 2D bipolarons will be strongly local- 
ized. Such 2D large bipolarons will unlikely be dissoci- 
ated in underdoped and optimally doped cuprates. Ac- 
cording to Emir* 5 - ' 39 ' 64 ' 65 , quasi-2D large bipolarons are 
attracted to one another via the phonon-mediated in- 
teraction that fosters their condensation into a normal 
liquid. Carriers are assumed to condense into droplets 
and form a liquid prior to their condensing further into 
a superconducting stated. The concentration of such 
large bipolarons above which they condense into a nor- 
mal liquid at a reasonable temperature is assumed to 
be xj;™' ~ O.l/nw^Ej where riu is the number of unit 
cells involved in a large bipolaron. To form one bipo- 
laron, two carriers are required from two dopants, then 
the critical doping level is x c = 2x£™ . Envisioning a 
large bipolaron involving about five cells, one cell and 
its four neighbors, the critical carrier concentration for 
forming a large-bipolaron liquid is about x c ~ 0.2/5 = 
q_q^5j39_ jj. wag a j g0 SU ggestec£i22 that in the cuprates, 
superconductivity occurs when the dopant density ex- 
ceeds the density of a large bipolaron liquid, namely, 
at x c > 0.05, and large-bipolaron superconductivity of 



the cuprates is determined by the competition between 
largc-bipolarons' attraction for one another and their at- 
traction to dopants. At the same time the condensa- 
tion of attracting large bipolarons (including also Cooper 
pairs of large polarons) into a superfluid Bose-liquid in 
the cuprates below the superconducting transition tem- 
perature T c was sugg ested in Refs^i 6 !. The conden- 
sation of large bipolarons into a superfluid Bose-liquid 
at sufficiently low temperatures in the cuprates can be 
regarded as an insulator-to-superconductor transition 28 . 
Emin argue d 5 ' 39 ! 64 ' 65 that (bi)polarons in the cuprates 
order in a manner commensurate with the underlying 
lattice, thereby forming a superlattice. For the square 
structure of CuC>2 planes such a commensurate ordering 
of large bipolarons can occur at x — 1/8 and lead to the 
conversion of a cuprate superconductor to an insulator 
composed of ordered large bipolarons 68 . Thus, the type 
of transition proposed by Emin is another possibility of 
the metal-to-insulator (at T > T c ) or superconductor-to- 
insulator (at T < T c ) transition which is driven by the 
commensurate ordering of quasi-2D large bipolarons and 
presumed superlattice formation. It is believed 39 that 
such transitions are neither a Mott transition nor an An- 
derson transition. 



V. CONCLUSIONS 

In this paper we have studied the carrier localiza- 
tion, MITs and stripe formation in inhomogeneous hole- 
doped cuprates driven by the gross electronic inhomo- 
geneities, charge ordering and extrinsic and/or intrinsic 
self-trapping of hole carriers. To model the real situa- 
tion in cuprates, the hole carriers are presumed to exist 
in a 3D ionic continuum medium and interact with the 
dopants (or impurities) via short and long-range carrier- 
defect-lattice interactions and with the lattice vibrations 
via short and long-range carrier-lattice interactions. We 
have investigated the self-trapping of carriers near the 
defects and away from the defects (i.e., in a defect- free 
deformable lattice) within the continuum model and adi- 
abatic approximation. Then, the possibility of the for- 
mation of extrinsic self-trapped states (i.e., hydrogenic 
and non-hydrogenic acceptor states) and intrinsic ones 
(large-polaronic states) in the CT gap of the cuprates is 
examined. We think that the large-radius dopants with 
the short-range repulsive potential can form hydrogenic 
impurity centers in LSCO and LBCO, where the hole 
carriers liberated from such impurity centers are self- 
trapped in a defect-free deformable lattice due to the 
strong short and long-range carrier-lattice interactions. 
When such hydrogenic impurity centers are distributed 
inhomogeneously and ordered differently with the for- 
mation of distinctly different superlattices, the narrow 
Hubbard impurity bands are formed in the CT gap of 
the cuprates. We have analyzed the applicability of the 
newly derived criteria for the Mott and Anderson MITs 
to the La-based cuprates. It is found that such MITs in 
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the cuprates with the hydrogenic impurity centers are un- 
likely possible. The situation in La-based cuprates with 
small-radius dopants is quite different from LSCO and 
LBCO cases. In these systems, the short- and long-range 
parts of the small-radius dopant potential are attractive 
and the combined effect of the strong carrier-defect and 
carrier-lattice interactions fosters carrier localization (or 
defect- assisted self-trapping of carriers) and formation of 
extrinsic large polarons. We have calculated the ground 
state energies (i.e., binding energies) of extrinsic and in- 
trinsic large polarons in La-based cuprates. The obtained 
results are consistent with experimental results on po- 
laroid formation in these systems. 

Further, we have demonstrated that the driving forces 
for carrier localization, MITs and stripe formation in 
La-based cuprates are the extrinsic and/or intrinsic self- 
trapping of hole carriers, the dopant and carrier-driven 
electronic inhomogeneities, the charge ordering and the 
anisotropy of dielectric constants. Under the certain con- 
ditions, the narrow energy bands of extrinsic and intrin- 
sic large polarons are assumed to prelude itinerant mo- 
tion of such polarons. These conditions and quantitative 
criteria for the new MITs in inhomogeneous La-based 
cuprates are obtained by using the pertinent uncertainty 
relation. We showed that the new MITs are accompa- 
nied by the formation of static and dynamic stripes, and 
would occur at different critical carrier concentrations 
depending on the ionic size of dopants (or impurities), 
the anisotropy of dielectric constants and the type of the 
ordering of extrinsic and/or intrinsic large polarons. In 
particular, the ordering of 3D intrinsic large polarons, 
the formation of different polaronic superlattices (z — 6 
and 8) and the anisotropy of dielectric constants favor 
the MITs and stripe formation in the range of doping 
x = 0.04— 0.125 in La-based cuprates with large-radius 
dopants. We argue that the x = 1/8 anomaly is charac- 
teristic of these cuprates, where the suppression of super- 
conductivity due to the formation of static stripes should 
be especially pronounced for the doping level x = 1/8. 
While the effects of the substitution of smaller cations 
for host lattice ions or dopants (e.g., Nd 3+ in LSCO), 
the anisotropy of dielectric constants and the ordering of 
extrinsic large polarons with the formation of different 
superlattices {z = 6 and 8) extend the metal-insulator 



crossover region from the heavily underdoped (x = 0.04) 
to the overdoped (x > 0.2) regime. Our results suggest 
that the new MITs and the formation of static stripes (in 
carrier-poor domains) and dynamic ones (in carrier-rich 
domains) occur in a wide range of doping (including also 
the "magic" doping 1/8) in La-based cuprates containing 
small-radius dopants (e.g., Ca 2+ and Nd 3+ ions), where 
the hole concentration x = 1/8 is not already particu- 
lar doping level corresponding to the formation of static 
stripes. 

When 3D extrinsic and intrinsic large bipolarons are 
formed in the cuprates at 77 < r/ c , such bipolarons will 
dissociate first into two large polarons at low carrier con- 
centration and then these polarons become mobile or ac- 
quire itineracy at high doping levels x > x c ~ 0.05. The 
essence of the above driving forces for carrier localization, 
MITs and static stripe formation in La-based cuprates 
has been mimicked by the theoretical predictions which 
capture the main characteristics of experiment. Thus, 
we have succeeded in explaining the experimental re- 
sults on carrier localization, MITs and stripe formation 
in hole-doped La-based cuprates with large- and small- 
radius dopants in the two characteristic doping regimes 
0.05 < x < 0.125 and 0.05 < x < 0.25, respectively. 

Finally, it should be noted that carrier localization and 
metal (or superconductor)-to-insulator transition accom- 
panied by the 1 /8 anomaly in hole-doped cuprates might 
be also driven by the commensurate ordering of quasi-2D 
(or planar) large bipolarons which will condense into a 
liquid above some critical carrier concentration. 
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